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Abstract 

We consider CP'^^^ models in d + 1 dimensions around Lifshitz fixed points with dynamical 
critical exponent z, in the large-N expansion. It is shown that these models are asymptocially 
free and dynamically generate a mass for the CP^~^ fields for all d = z. We demonstrate 
that, for z = d = 2, the initially nondynamical gauge field acquires kinetic terms in a way 
similar to usual CP^~^ models in 1 + 1 dimensions. Lorentz invariance emerges generically in 
the low-energy electrodynamics, with a nontrivial dielectric constant given by the inverse mass 
gap and a magnetic permeability which has a logarithmic dependence on scale. At a special 
multicritical point, the low-energy electrodynamics also has z = 2, and an essentially singular 
dependence of the effective action on B = eijdiAj. 



Nonlinear sigma models are ubiquitous in a variety of areas in theoretical physics. In this 
paper we will deal with the CP^'^ model [1], whose fields are N component complex vectors 

— * 

x) constrained by 

(1) 
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and fields which differ by an overall (space-time dependent) phase are identified, 

${t, x) ~ e^^(*'^) x) (2) 

The identification is incorporated by introducing a non-dynamical U{1) gauge field, A^. The 
conventional, relativistic, action is 

S=]^j dtj (fx {Djy{D>^$) (3) 

where 

D^ = d^ + lA^ (4) 

— * 

Integrating out A^t, leads to a nonlinear action which involves only the fields. As is well 
known, m. d — \ the model (3) is asymptotically free and generates a mass m for the fields 

— * 

by dimensional transmutation, as can be explicitly seen in the 't Hooft large-N expansion 
N ^ oo g ^ g^N — X — fixed [1]. At the same time, the initially nondynamical 
gauge field acquires a standard kinetic energy term, with a gauge coupling constant given by 
m^. This is the simplest example of a dynamical emergence of gauge dynamics. 

The d — 2 model is interesting for condensed matter applications. In fact, the 0(3) nonhnear 
sigma model with three component unit vector n can be rewritten as the CP^ model via the 
identification 0(T0 = n. It is evident that local phase transformations of the fields do not 
affect the "gauge-invariant" field n. Now there is a usual order-disorder transition: In the 
magnetically ordered phase of n, the field is condensed, and gauge field is gapped out by the 
Higgs mechanism. However, gauge field dynamics appears in the disordered phase when the 
(j) fields become massive. Normally, the gauge fields also become massive and fields become 
confined on the paramagnetic side of the transition due to the compactness of the gauge field 
[2]. However, suppressing the monopoles [3] of the gauge field (which correspond to "hedgehog" 
configurations of the original n fields) leads to a new critical point [4] , and a paramagnetic phase 
with a gapless photon [5] . There are conjectures that such a model with a noncompact gauge 
field also describes a possible non-Landau, deconfined, critical point [6] between the Neel and 
bond-ordered phases of the d — 2 quantum antiferromagnet. 

In this paper we consider UV modifications of these models, which correspond to Lifshitz-like 
fixed points with a dynamical critical exponent z, 

SL^^JdtJ d^x [(L>o0)*(dV) + a(A0)*(^V) + l^'^r] (5) 



where the operator T>^ is a sum of 0{d) invariant terms containing z factors of the spatial 
covariant derivative D^. For example, in 2; = 2 



V^ct>\ 



2 _ 



(6) 



with a, 6 > being parameters. For higher z we would have many more terms corresponding 
to various orderings of the A- 

At the fixed point a = 0, one needs to scale the time and space coordinates as 



Such fixed points, called Lifshitz fixed points, have a variety of applications in classical con- 
densed matter systems [7]. They also have a connection to quantum dimer models [8], which 
are defined with a "kinetic" term which flips dimers on parallel bonds, and a "potential" term 
which gives an energy to every flippable plaquette. Finally, there is a constraint that every lat- 
tice site should have one and only one dimer touching it. Recently, it was realized [9, 10] that 
neutral (ungaugcd) onc-componcnt Lifshitz fixed points describe special points (the Rokhsar- 
Kivelson (RK) points) of quantum dimer models [8] on bipartite lattices, where the field is a 
height variable dual to a bond which may or may not contain a dimer [11, 12]. The standard 
RK point describes the transition between the smooth and rough phases of the height and is 
multicritical, in the sense that more than one parameter needs to be tuned to attain the fixed 
point [9, 10]. However, it is possible to construct models with enough symmetries such that 
the fixed point can be obtained as a regular crtical point describing, for example, the phase 
transition between two different types of bond-ordered states in a bilayer honeycomb latticc[13]. 
While the above are examples ol z = d = 2 theories of neutral scalars, examples oi z = 2 gauge 
theories also occur in condensed matter, in the description of algebraic spin liquids in 0? = 3 
[14, 15] and topological critical phases in d = 2[16]. 

Note that the action of Eq. (5) still contains a single time-derivative of the gauge fields Aj, 
though it has higher spatial derivatives. Thus, even though one cannot easily integrate out the 
Afj^ to obtain a pure spin model, the gauge field is non-dynamical to begin with. In the spirit 
of the renormalization group, one expects that the model defined by Eq. (5) for A'" = 2 is in 
the universality class of a z = c? = 2 0(3) nonlinear sigma model. 

Recently, Lifshitz-type theories have been suggested as UV completions of low energy 
Lorentz invariant theories of gravity and gauge dynamics [17]. This is because theories which 
are non-renormalizable at the usual Lorentz-invariant UV fixed points with z = 1 can become 
renormalizable for non-trivial z. The idea that Lorentz symmetry violation can be regarded 
as UV regulators of field theories has been around for a while. See [18] for a recent discus- 
sion and references. In the same spirit, recently such Lifshitz fixed points have been proposed 



t 



X — > 7i 



(7) 
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as UV completions of four fermion theories similar to the Nambu-Jona Lasinio model in 3+1 
dimensions and discussed as possible candidates for physics at the weak scale [19], [20]. 

We will find that for any choice of the theory defined by (6) is asymptotically free for 
all d — z and generates a mass gap, pretty much as the models in [20] - so that the theory is 
always in a disordered phase. As a result, the gauge field acquire a kinetic term. 

This means that the possible emergence of Lorentz invariance at low energies is a little more 
non-trivial since this has to happen in the gauge as well in the scalar sector. For z — d — 2 
we explicitly show that this indeed happens generically, and discuss the possibilities for higher 
d — z. 

For special choices of T>^ (the a = multicritical point in our case) something even more 
interesting happens: one obtains 8k z — 2 electrodynamics with a standard E'^ term, but the 
leading term m B — eijdiAj is (V-B)^. This naively suggests that a constant B costs no 
energy: however a more careful calculation reveals that there is a nonanalytic dependence on 
constant B of the form S^/^ exp {—7rm/B). Note that the analytic terms in our z — 2, d — 2 
electrodynamics have the same form as the gauge theory descriptions of algebraic spin liquids 
in d — 3 [14, 16], but appear to be dual to the gauge description of the transition between two 
bond-ordered phases [13] or the topological critical phase[16] in d = 2 in which it is the E'^ 
term which is replaced by {eijdiEjY. 

1 Asymptotic Freedom and Dynamical Mass Generation 

We will study the large-N hmit of the model (5) with the constraint (1), using standard tech- 
niques [21]. The coupling g in the model (5) becomes dimensionless under Lifshitz scahng (7) 
when d — z. This may be seen by dimension counting : t has length dimension z, so that the 

— * 

length dimensions of and g are 



Whether the coupling is marginally relevant or marginally irrelevant ai z — d depends on the 
dynamics. To investigate this we use standard large-N techniques. Imposing the constraint (1) 
by a largrange multiplier field x(i, x) we get the action 



[0]^[L]^ 



[g] - [L]'2 



(8) 



^L^lJdtJ d^x {Do$f{D% + a{D,$Y{D'$) + \V'$f + x(t, x) [\$\^{t, x) 



1 



) 



(9) 



— * 

Integrating out the field we get the effective action 



Seff = 7V{Tr log [-Dl + {-iy{V^f + x(t, ^)] - ^ / dtd^x x{t, x)} 



(10) 
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At A?" = oo the functional integral over Ai^{t,x) and x(^)^) is dominated by the saddle point 
of (10). We will assume that the saddle point is translationally invariant and rotationally 
symmetric in the d spatial dimensions with a vanishing gauge field strength. Thus we may set 
x(i,x) — Xo in the saddle point equation ^. This also means that so far as the saddle point 
equation is concerned, all possible terms in contribute equally 

2N I I = 1 (11) 

For any finite a, this integral is logarithmically divergent ior d — z and behaves as log^, 

— * 

where A is a cutoff on the spatial momentum k. This immediately implies that a solution to 
the gap equation is 

m2 = Xo~A2^exp [--^] (12) 

where A is a positive real number . Since m is (to leading order in 1/A^) the physical mass 
of the (j) field (i.e. in a lorentzian signature this is the lowest value of the energy of a single 
particle state), it is clear from (12) that the coupling has to be asymptotically free, with a 
beta function 

* d g^N 

It is useful to evaluate the integral in (11) for our primary case of interest, z — d — 2. For 
a, m <^ A^ we get 

n —-52- Ct 

m = 2A^ e ^ - ^ (14) 

The standard gaussian fixed point corresponds to a » A^ and leads to a linearly divergent 
answer in this case. 

Dynamical mass generation for this model is thus almost exactly identical to that in the 
four-fermion model of Ref . [20] . The effective action for the gauge field and the fiuctuations 
of x(i, x) has to be now obatined by substituting 

X{t,x) ^Xo + A^{t,x) ^ -^A^{t,x) (15) 



in (10). Clearly, this will generate kinetic terms for Sx and A^. The effect of these will be to 
provide corrections to the leading order propagator of the fields which is simply the integrand 
of (11). Accordingly, the parameter a will be renormalized. If we go off the critical surface 
containing the Lifshitz fixed point, the renormalized value of a will be nonzero. Clearly, when 



^In principle there could a condensation of the field strength. However we will soon see in Section 2.1 that 
ioi d = z = 2 the effective action for a constant B = CijdiAj is always larger than the action with B = 0. 
This rules out condensation of B. Our assumption that the field strength vanishes at the saddle point is thus 
justified only a posteriori. We do not have a proof that this continues to hold for all d = z, but this appears to 
be plausible. 
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— * 

the spatial momenta are much smaller than -^/a, the propagator of the (f) will be dominated by 
the arenk'^ term. Therefore at low energies when a 0, Lorentz invariance is recovered with a 
speed of light given by 

2 Effective Action for the Gauge Fields : d = z = 2, a = 

Emergence of Lorentz symmetry at low energies in the gauge field sector is more non-trivial, 
especially when a — 0, which is the case we will concentrate on. By gauge invariance, the 
induced action for the gauge fields must be functionals of the field strengths Foi and Fij and 
their derivatives. In addition, it must be symmetric under spatial rotations. For a Lorentz 
symmetry to emerge, this effective action must contain combinations like 

^oF^ + (16) 

with constant eo, Ho- In that one can now rescale t, x, Aq, Ai to get a standard Lorentz invariant 
form 

The length dimensions of the dielectric constant eg and magnetic permeability /Iq may be 
easily seen to be 

[6o] ~ [Lf-"^' [f^o] - [L]'^'-' (17) 
so that the speed of light c — 1/ yZ/Io^ has length dimensions 

[c] ~ [L]'-^ (18) 

as it should. 

It is not at all obvious that terms like (16) have to emerge at a = 0, since the parent theory 
has z — 2. In fact we will show that for special choices of the operator T>^ this will not happen. 
However, we will find that for generic choices of T>^, terms hke (16) do appear. 

Let us first address this question ior z — d — 2, using the form (6). For this purpose, it is 
sufficient to consider the effective action (10) with 5x — 0, so that we essentially have 

Seff = N{Ti log [-Dl + {a + h){DiD'f + a{B{t,x)f - iae'^{diB)Dj + m^]-{B = term)} 

(19) 

where we have used the commutation relation 

[A, Dj] = iFij (20) 

and for d = 2 renamed F12 — B{t, x). 
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2.1 Constant Magnetic Field 

It is useful to first evaluate this for a constant B. Then the problem in evaluating the effective 
action reduces to the problem of determining the eigenvalues of the operator 

H(B) = -Dl + (a + h){-Dl - Dlf + aS^ (21) 

which is closely related to the problem of Landau diamagnetism. Let us choose a Landau gauge 

Aq^A^^Q A2^B (22) 

Consider the system to be in a large box with size in the time direction T and spatial sizes 
Li,L2. For large enough T,Li,L2 the eigenvalue of do can be taken to be continuous, which 
we will call Pq. It is straightforward to see that the eigenvalues of H{B) are 

«(Po, n) = + (a + h)B'^{2n + 1)^ + aS^ (23) 
with a degeneracy of the level n given by 

d(n) = ^ (24) 
To evaluate the effective action (19) we use the Nambu-Schwinger-de Witt representation, 

-Seff ^nJ^^ e-"""' TV e-^^(^) - (S = term) (25) 
Using (24) and (23) we have 

n ^-00 2tt 2tt 

e-^«^' i?2 [0 I 4iB^s{a + b)/7r] 



n=0 

VT B fW 



87r2 



(26) 



where V — L1L2 denotes the spatial volume and '(92[w|r] is a Jacobi thcta function To 
examine the small B behavior it is useful to use standard theta function identities to write 

VT 

rj^ ^-sH{B) ^ ^-saB^ I ,7^/(45^5(0 + h))] 

IGtts \/a-\-h 

_..R2 ^ ^_^^k g-4.B^(„+,) (27) 
fe=— 00 





IGtts ya + h 



The theta function in (27) may be written in a product represenation as 

i?4 [0 I ml{AB^s{a + h))] = - e «^ 2 (1 _ g iB^s(a+b) ) (28) 



n=l 



^We are grateful to Al Shapere for pointing out that an efficient way to manipulate this sum is to recognize 
this as a theta function. 
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Since a,b >0, this immediately shows that 



Tr e-^^(^) < Tr e'^^^") (29) 

so that 

SeffiB) > S^ffiO) (30) 

for any m^. This provides a justification for setting = in the saddle point equation which 
determines m^. The result (30) in fact holds for all d with z = 2. 

The integral over s in (27) can be performed, leading to the effective action ^ 



SeffjB) ^ ^Bm aB-^ TTkm 



1 + ^^ 


1 


a + b 





(31) 

where Ki denotes a modified Bessel function. In deriving (31) we have noted that the k — 
term in the sum in (27) is the sole contribution when B — and subtracted that. 

The first term on the right hand side of (31) has a non-analytic dependence on B for small 
B. This follows from the asymptotic behavior of the modified Bessel function. For B <^m the 
sum in (31) is dominated by the k = 1 term, which leads to 

SeffjB) - SefM _ B-2 vnh b\ 

The second term contains various powers of B. 

Therefore, at the multicritical point a = 0, the effective action vanishes for S = in a 
non- analytic fashion ^. 

For any a 7^ the final form of the effective action begins with a term proportional to S^, 
It follows from (25) that the coefficient of this term is divergent (proportional to r(0)). To 
understand this, we use a version of dimensional regularization by adding {d — 2) extra spatial 
directions. Now the first line of (26) will be modified to 

jj_Q J— 00 27r J (27r)'^ ^ 27r 



This leads to the following coefficient of B^ in the effective action (25) 

ds d-2 r(^i 



„r^£e-»'.,-^=„ii5J (34) 

Jo S 777, 2 



^Note that there is an overall factor of A'' in the effective action (39). Since we are performing a 1/A'' 
expansion, so that the fields have to be rescaled as in (15). The factor of N cancels for the terms which are 
quadratic in the fields. 

^We have performed the sum in (27) using a Euler-McLaurin expansion and verified that there are no 
polynomial terms in B in this case to very high orders. 



The coefficient lias lengtfi dimension (2 — d) as required. In a dimensional regularization sdieme 
we introduce a scale k with length dimensions —2, and write this coefficient as 

1 



2-d 



(35) 



where is dimensionless. Then as e — > 0, the finite part of 1/no becomes 



^ log (-) (36) 



Ho \mj 

2.2 General Calculation Using Heat Kernel 

In this section we will evaluate the effective action using a standard heat kernel method for 
general fields -Foi(t, x) and Fij{t, x). The calculation presented above for constant magnetic field 
shows that at the special point a = the effective action has a non-analytic dependence on B. 
We would like to determine whether the action contains terms which are analytic in derivatives 
oi B. As shown in the Appendix , the small s expansion of the heat kernel is of the form 

/oo 
d^xdtY,hn{x) (37) 
n=l 

which leads to the effective action 

Seff ^NJ2 d^^dt hn{x) (38) 

As in the previous subsection the term with n = 4 can be handled via dimensional regulariza- 
tion. As explained in the appendix, we should therefore replace (38) by 

oo . Yi- - 1 -) 

Seff ^NY, I d^xdt bn{x) ^'2(^-1+1) (39) 



n=l 



Let us first evaluate this for a = and 6=1 For this case, explicit calculations yield bi{t,x) = 
b2{t,x) = bs{t,x) = b4{t,x) = b^{t,x) = 0. The leading contribution comes from bQ{t,x). After 
a fairly long calculation we find that this leads to the effective action 

1 



S,ff,a^O = ^Jdtd^X 



Fl + Y^^d^BWB) 



(40) 



Higher powers of field strength are suppressed by powers of N. Note that the term we ob- 
tained in the previous section for constant B is nonanalytic in B, and formally irrelevant by 
power counting since its Taylor series is 0. However, in the absence of this nonanalytic term, 
any constant B costs no energy, leading to a huge ground state degeneracy. Thus, the term 
^3/2g-7rm/B jg ^ daugcrously irrelevant operator for this special case a = 0, 6 = 1. Clearly 
Lorentz invariance is not regained at low energies in this case. 
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Away from this multicritical point a 7^ 0. As expected from our constant B calculation we 
now find that 64 7^ 0. It is straightforward to see that 

h{t,x)^r^Fij (41) 

Since we are interested in the action at low energies, we should retain only the lowest non-trivial 
terms with the least number of derivatives. This leads to the low energy effective action for 
gauge fields, upto numerical factors 

Seff ~ / d'^x dt [-Fo,F°^ + ^FijF'^ + ■■■] (42) 

where the ellipses now stand for terms containing more derivatives and/or more powers of the 
field strength. As in the previous subsection, in the spirit of dimensional regularization this is 
really 

S,ff ~ J d^x dt [^^oii^°' + log F,^F'^ + • • •] (43) 
so that we have /xq given in (35) and 

eo (44) 

m 

Since eo, //o are constants independent of (i, x) one can now rescale x, Aq, Ai to get the form 

5e// ~ j d^x dt + ■ • •] (45) 

This demonstrates the emergence of approximate low energy Lorentz symmetry with a scale 
dependent speed of light. Note that this happens even when the (renormalized) parameter aren 
in (5) is tuned to zero. At this Lifshitz point there is no lorentz symmetry in the scalar sector. 
For a^en 7^ the speed of light in the scalar sector is different from that in the gauge sector. 

3 Discussion 

It is clear from section (1) that the coupling g in the sigma model is asymtptically free for all 
d — z. Does this also mean that gauge dynamics emerges in all dimensions? Prom equation 
(17) we see that the length dimensions of eo is always 2 for all z — d. If such a term appears, 
one would expect that eo ~ m^/^, since from (12) the length dimension of the mass m is 2;. For 
2; > 3 it is not clear how such a term can arise in the effective action obtained by integrating 
out the massive field 0. In fact if such an effective action is an expansion of powers of 1/m 
(apart from logs) one would expect that the lowest dimension operator which would appear 
must have length dimensions [L]^''"''. In a similar vein, the magnetic permeability would have 
positive length dimensions. Such a term is also unlikely to come from an effective action. It 
would be interesting to investigate this issue further. 
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The model studied in this paper gauges the overall U (1) of the symmetry group. One could, 
instead, consider gauging the entire U (N) group to obtain a non-abelian gauged sigma model. 
This model would generate a mass gap in exactly the same fashion - in fact the gap equation is 
identical. It would be interesting to see the effective action for the non-abelian gauge fields in 
this case. It appears to us that the heat kernel expansion calculation is quite similar to ours. 

Another interesting direction is to revisit the linear rather than the non-linear sigma model 
around a Lifshitz fixed point as has been originally considered in classical statistical mechanics. 
The length dimension of a (0 • 0)^ couphng is given by {d — 3z), so that this is a relevant 
operator for z > d/3. It would be interesting to explore if there are IR fixed points for d > 2 at 
finite values of the corresponding coupling, similar to z — l,d — 2. This could have interesting 
applications to particle physics. These vector models can be also interesting from the point of 
view of AdS/CFT correspondence. Lifshitz fixed points have been argued to have dual gravity 
descriptions [23]. On the other hand, the dual of usual vector models are higher spin gauge 
fields in usual AdS [24], [25]. It would be interesting to see the nature of the gravity duals for 
these Lifshitz sigma models. 

These issues are currently under investigation. 
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A Expansion of the Heat Kernel 

In this appendix calculate the heat kernel using the technique of [22]. This uses the represen- 
tation 




(46) 



Using the basic identity 



(47) 



we have for our case 



^ u;^ + + ai + a2 + a3 + a4 



(48) 



where 
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as = i[{D'Di){kWj) + {kWj){D'Di) + 2Di{kWj)D'] 

04 = -Dl + DiDjD'D^ (49) 



Using (48) in (46) and rescaling 



1 - 1 - 



we get 
where 

113 

G = s^a-^ + 52^2 + s^Qg + s a4 (52) 

— * 

The integrals over u and k can be now evaluated, leading to small-s expansion of the heat 
kernel of the form (37), leading to the form of the effective action (38). 

The term with n = 4 has to be treated in dimensional regularization. This means that in 
(46) we replace ^k — > df'k^ so that after the rescalings (50), the equation (51) becomes 

r^^-sHiE,B) ^ _L / dtd^^ \ e-(-'+^') e-«(^'^) (53) 

s— J J (2tiV 



For d = 2 — e we will still evaluate the integrals over k by replacing d'^k — > d'^k in (53). This 
leads to the small s expansion 

/oo 
d^xdt^bnix) st-i+t (54) 
n=l 

This leads to the expression (39). 

The integrals over lu and k can be performed using basic symmetry properties. Thus, e.g. 

/oo „ 

dcu e-'^ 0;^"+^ = (n integer) 
-oo 

/oo 2 o 1 2 

duj uj —7: 
-00 Z J —00 

/oo 2 3 2 

duj e~'^ uj"^ — - j duj e""^ 
-00 4 J— 00 



(55) 



while 



J d^k k'^k'^ k^^-+' e-''' = (n integer) 

/ 



d^k k'y 6-^=' = -^6'^ / d'k e 
2^ 
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/ 



j£k{k-kf e-^" = l^Jd'^ke-'^' 

cPk ik ■ kf k'y e-^" = -^5'^ I cPk e-^' 

2y TT J 



(56) 



Using these integrals, it is straightforward to see that terms with n = 1, 3, 5 in the sum (54) 
vanish since they have odd numbers of u and/or k. A short calculation using the explicit 
expressions in (55) and (56) shows that the various terms cancel, leading to b2{x,t) = 0. The 
first non-trivial term is therefore for n = 4. Here, after several cancellations one is left with 

h(t,x)r^FijF'^ (57) 

which basically comes from rewriting the second term in 04 as 

DiDjD'D^ = DiD^Di+DiDj[D\D^] = DiD^Di + ^[Di, Dj][D\ D^] = DiD^Di + ^F^jF'^ (58) 

The next nonzero term comes at n = 6. This leads to an electric field term, which arises from 



dujdfik 



(01020102) (59) 



(27r)3 

which clearly includes a term 

DiDoDiDo (60) 

and hence to FqiF^. The n — 6 term contains other contributions as well. These contain higher 
derivative terms in the field strength B. Specifically, for the case a — 0, b — 1, we obtain the 
term 

-^(VBf (61) 
120m^ ' ^ ' 
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